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Abstract 

We present the results for a perturbative determination of mass dependent improvement coeffi- 
cients is, r s , ce and cb in a relativistic heavy quark action, which we have designed to control niQa 
errors by extending the on-shell 0(a) improvement program to the case of uiq 3> Aqcd> where 
niQ is the heavy quark mass. The parameters v and r s are determined from the quark propa- 
gator and ce and cb are from the on-shell quark-quark scattering amplitude. We show that all 
the parameters, together with the quark wave function and the mass renormalization factors, are 
determined free from infrared divergences once their tree level values are properly tuned. Results 
of these parameters are shown as a function of niQa for various improved gauge actions. 
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I. INTRODUCTION 

A calculation of weak matrix elements for the B and D mesons is a subject of great 
interest in lattice QCD: their precise determination is an essential ingredient to extract 
the Cabibbo-Kobayashi-Maskawa matrix. Although in principle we can extract these weak 
matrix elements precisely from lattice QCD simulations, it is still difficult to achieve this 
goal. A main obstacle is the systematic error originating from large mga corrections with 
current accessible computational resources: m^a ~ 1 — 2 and m c a ~ 0.3 — 0.6 in the quenched 
approximation and m^a ~ 2 — 3 and m c a ~ 0.6 — 0.9 in unquenched QCD. 

To control this large mga corrections a new relativistic approach is proposed from the 

n 

view point of the on-shell 0(a) improvement program 1]. The generic quark action is given 



where we are allowed to choose r t = 1 and other four parameters v, r s , ce and cb are analytic 
functions of mga and the gauge coupling constant g. 

In this formulation the leading cutoff effects of order (mga) n are absorbed in the definition 
of renormalization factors for the quark mass and the wave function. After removing the 
next-leading cutoff effects of 0((mga) n a Aqcd) with u, r s , ce and cb in the quark action 
properly adjusted in the rriQa dependent way, we are left with at most 0((aAQCT)) 2 ) errors. 

In this paper we calculate u, r s , ce and eg up to the one-loop level for various improved 
gauge actions. While v and r s , together with the quark wave function renormalization 
factor Z q and the pole mass m p as byproducts, are obtained from the quark propagator, 
ce and cb are determined from the on-shell quark-quark scattering amplitude. One- loop 
diagrams are evaluated by employing the conventional perturbative method with the use 
of the fictitious gluon mass to regularize the infrared divergence, which was successfully 
applied in the massless case to the calculation of the renormalization constants and the 
improvement coefficients for the quark bilinear operators Q| and the improvement coefficient 
cswjsl- We show that the improvement parameters u, r s , ce and cb are determined free 



by 



S q = m q(x)q(x) + q(x)-f D q(x) + vJ2q(x)-fiDiq(x) 





(1) 
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from the infrared divergences, once their tree level values are correctly tuned in the rriQa 
dependent way. Some preliminary results are presented in Ref. jj]. 

This paper is organized as follows. In Sec. II we introduce the relativistic heavy quark 
action and the improved gauge actions. We also give their Feynman rules relevant for the 
present calculation. We describe the method for the calculation of the one-loop diagrams in 
Sec. III. In Sec. IV we determine v and r s from the quark propagator. The results for Z q 
and m v are also presented. In Sec. V c# and ce are determined from the on-shell quark- 
quark scattering amplitude. We give a rather detailed description on the cancellation of 
the infrared divergences in one-loop diagrams. In Sec. VI we explain how to implement the 
mean field improvement of the parameters. Our conclusions are summarized in Sec. VII. 

The physical quantities are expressed in lattice units and the lattice spacing a is sup- 
pressed unless necessary. We take SU(iV c ) gauge group with the gauge coupling constant 
9- 



II. ACTION AND FEYNMAN RULES 

The relativistic heavy quark action proposed in Ref. jl| is given by 

Squark = ^ ~ {V>»(-?"t + 7o) ^,oVVf6 + ^H"* ~ 7o) ^-^(A-o} 
n 

+ 12 o 12 {^n(-r s + vji)U n ^ n+i + M-r s - viiW^-. .ifj n _-} 

n Z i 

+(m + r t + 3r s ) ^ ip n ip n 

n 

-^1212^ '^n^OiFu{n)^ n - CB^^^J^nCTijFij^lpn, (2) 
n i ^ n i.j 4 

where we define the Euclidean gamma matrices in terms of the Minkowski ones in the 
Bjorken-Drell convention: 70 = j BD , 7, = -rf BD {j = 1, 2, 3), 75 = ^ BD and = §[7 M , 7J. 
Whereas the value of r t can be chosen arbitrarily, z/, r s , c B and c B have to be adjusted 
to remove the cutoff effects of O ((mga) n a Aq C d)- As explained in Ref. jl| the (mga)" 
corrections can be avoided by the redefinition of the quark field and mass. The field strength 
F^y in the clover term is expressed as 

F,M = l£j-ton)-0?(n)), (3) 
4 i=i M 9 

U x (n) = Un^Un+^Ul+^Ul,, (4) 
3 



U 2 (n) - UnpU^+fi^Ul^Un-fcpt, 



n— fj,— w n— v,v i 



Un 



(5) 
(6) 
(7) 



The weak coupling perturbation theory is developed by writing the link variable in terms of 
the gauge potential 



U n!fl = exp [igaT A A A f n + 
N% — 1) is a generator of color SU(iV c 



where T A (A 



(8) 



The quark propagator is obtained by inverting the Wilson-Dirac operator in eq.©, 
Sq 1 ^) = «7osin(p ) + z/i^7isin(pi) + m 

i 

+r t (l - cos(po)) + r s ^(1 - cosQo;)), (9) 



For the present calculation, we need one-, two- and three-gluon vertices with quarks: 



Vw(p,Q) = -#T A jz7 cos 
V x A {p,q) = 



V** C ( P ,q) 



+ r t sm I — - — 
Pi + q% 



rrA J • I Pi + 1i i 

< ^7iCOS ( — - — ) + r s sm 



2 9 2 {T ' T >h oSm (-2--J- r « OOB (-2-- 

a 2 l i r r Ar r B X { ■ ■ (Pi + qi\ fPi + qi 

-g -{T ,T }{ vl7tBai \——)-r.ca B \ c - r - 
a\ 3 l 



—g~- T {T , T } + T {T , T } + T {T , T } 
6 6 L 

x jry cos I — - — j +r t sin I — - — 

a 2 J 
— o - 
6 y 6 

+ . • [Pi + qi 

x < wji cos I — - — ) + r s sm 



^(P.?) = -^-^c S w^cos 



/V - q^ 



sm{p v - q h 



9 2 if A BcT c 



X \ Csw^V^v 



1 

4 

/ fc^N / fc 2A A (q»-Pn , 

4 cos — cos — - cos — cos 

V 2 J \ 2 \ 2 



(10) 

(11) 
(12) 

(13) 

(14) 

(15) 
(16) 



qv -Pv 



-2 cos cos(^ 



(17) 



+<V csw^^p sin 



q^-Pi, 



[sin(fc 2p ) - sin(A; lp )] 



X 



T A T B T%, V J] zc sw MP ^ p {-i cos (^y^) sin(g p - p„) 



,'%-Pn\ (Q P -P P \ (hp- hp] . (hp 
( " s I — 2 — J COS V — 2 — ) C ° S I 2 / Sm I 2 



' - "-' icWV,,,, (18) 
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rpArj-iBrj-iC _j_ rpCrpBrpA 



x {^2 cob cos cos (^J^j sill 

+<W sm I I cos I — ^ + fc 2 ^ I 

. /fc lAl + 2k 2fl + k 3 A fq v -p v \ (hv-hu\ 
+ V sm ^ j cos (^^— j cos 1^^— J 

where Csw°* — c sw*° — c e, c sw 4,7 — c .b (z, j — 1)2,3) and /abc is the structure constant of 
SU(iV c ) gauge group. The first six vertices originate from the Wilson quark action and the 
last three from the clover term. The momentum assignments for the vertices are depicted 
in Fig. [TJ 

For the gauge action we consider the following general form including the standard pla- 
quette term and six-link loop terms: 

S S = \ | C H tlU pl + Ci tTU rtg + tiUchr + C 3 ^ txU p l g > (19) 

9 ^ plaquette rectangle chair parallelogram J 

with the normalization condition 

c + 8 Cl + 16c 2 + 8c 3 = 1, (20) 

where six-link loops are composed of a 1 x 2 rectangle, a bent 1x2 rectangle (chair) 
and a three-dimensional parallelogram. In this paper we consider the following choices: 
ci = c 2 = c 3 = O(Plaquette), C\ = —1/12, 02 = 03 = 0(Symanzik)|j], 0] ci = —0.331, 
c 2 = c 3 = O(Iwasaki), ci = -0.27, c 2 + c 3 = -0.04(Iwasaki') 0, c a = -0.252, c 2 + c 3 = 
-0.17(Wilson)3 and ci = -1.40686, c 2 = c 3 = 0(doubly blocked Wilson 2 (DBW2))fl. 
The last four cases are called the RG improved gauge action whose parameters are chosen 
to be the values suggested by approximate renormalization group analyses. Some of these 
actions are now getting widely used, since they realize continuum-like gauge field fluctuations 
better than the naive plaquette action at the same lattice spacing. 
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The free gluon propagator is derived in Ref. [5|: 

1 



with 



The matrix A M „ satisfies 



2\2 



1 ~~ Ap U )kpk u + 5p U 2_j k<,A v 



(i) = for all fi, 
(h) = A Ufl , 

(iii) A^{k) = A^i-k). 

(iv) A^{0) = 1 for n ^ v, 

and its expression is given by 

(k — k^)(q^ p q^ T kj l + q^pqprkp + qp T qprk T 
-\-{k — k p)(c[ v p q VT k v + q u pqprkp + qm-Qprkr) 
+q w q,r(kl + k 2 p ){k 2 v + k\) + + k 2 T )(k\ 



i Vp ' •"pJVV ' T/ 1 1P> I^PV" p ' 

(fc p + k T ) - {q^qup + qp T qur)k p k T 

ill 2 , ?2?2 , ?2?2 , i2 



q.pvQ.pTV p 

~Qpu(qpp K p^ T 



p t + Qprkpkp + q V pk v k T + qvrkjk^] 



with n ^ u ^ p ^ t the Lorentz indices, qy^ and A 4 are written as 



1 - 5p U ) [l - (ci - c 2 - c 3 )(fcj + fc*) _ ( C2 + C3 )^ 

"2r,2. 



A 4 = Z fc J II ?*7i + S ^pklqpv{q w qur + qprqup)- 

In the case of the standard plaquette action, the matrix v4 Mi , is simplified as 



/iplaquette i r 



The present calculation requires only the three-point vertex which is given in 



-i^fABC C i V 9 %r( k U fa, fa) 
8=0 
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with 

vff XpT (k 1: k 2 ,h) 



8\ P (h - k 2 ) r c 3X + 2 cycl. perms., 



(32) 



V$ XpT (h,k 2 ,k 3 ) = 16V£{ pT (h,k 2 ,k, 



+ 5 Xp {c 3A ((A;i - k 2 ) x (5x T kl - k 3X k 3r ) - (h - k 2 ) T {kl T + kl T )) 
+ (fci - k 2 ) T (k lx k 2X - 2c lx c 2X k 2 x )} + 2 cycl. perms. 

(0) 



(33) 



6x P (l-5 Xr )c 3X £ {(*i-W& + & + &) + feh-(&-&)} 



+ (1 - 5 Xp )(l - 5 Xt )(1 - 6p r )k 1 \k 2 p(ki - k 2 ) T + 2 cycl. perms. 

r(0) 

g3Ap? 



(34) 



SV^Lih, k 2 ,k. 



Sxp(l ~ hr)c 3X (h - k 2 ) T £ Ckl + kl) 
+ 2-(l - VX 1 - 5at)(1 - 5 pT ){k 1 ^k 2 ) T \k lx k 2p - \{k^k l ) p {k 2 ^k 3 ) x 



+2 cycl. perms.] 



where we introduce the notation, 



(35) 



CiX 



cos 



VA 



(36) 



The momentum assignments are found in Fig. El 



III. METHOD FOR THE CALCULATION OF ONE-LOOP DIAGRAMS 



We determine Z q = Zf> + g 2 Z^ 



q ) m P 



,2„(1) 



(0) + n 2 u (l) 



M 1 ) 



from the quark propagator and eg = c§ + g 2 c E l , cb = Cg + g 2 Cg from the on-shell 
quark-quark scattering amplitude, where the superscript (i) denotes the i-th loop level. 
The tree level values of the parameters are already determined in Ref. To calculate 
the one-loop contribution we write a Mathematica program for a given loop diagram to 
compose the integrand of the Feynman rules and to practice the Dirac algebra. The output 
is then transformed into a FORTRAN code by Mathematica. The momentum integration 
is performed by a mode sum for a periodic box of a size L 4 with L = 64 after transforming 
the momentum variable through k'^ = k^ — sin/c M . The numerical errors are estimated by 
varying L from 64 to 60. 



Although z/ 1 ), r^-\ c# and Cg should be determined free from the infrared divergence, 
we have to deal with it in the process of calculations. Following the method employed 
in Ref. ^(|, we subtract, from the original lattice integrand of each one-loop diagram, a 
continuum-like integrand which is an analytically integrable expression with the same in- 
frared behavior of the lattice integrand. Suppose -4^4 1 (k, {p}, A) is a lattice Feynman 
integral for a given one-loop diagram, whose infrared divergence is regularized by the ficti- 
tious gluon mass A. {p} denotes a set of external momenta. The infrared divergent term is 
extracted as 



A->0 



+ Lwf i(kM ' (37) 

where I(k, {p}, A), which has the same infrared behavior with I(k, {p}, A), should be ana- 
lytically integrable. The infrared divergence is transferred in the last term. A candidate of 
the counter integrand I(k, {p}, A) depends on each one-loop diagram. We will explain it in 
the following sections. 



IV. DETERMINATION OF v AND r, AT THE ONE-LOOP LEVEL 



At the tree level the parameters are adjusted such that the quark propagator of eq.(jHJ) 
reproduces the correct relativistic formQ]: 



1 -nWo - i Ei liPi + 2 

- + (no pole terms) + (J{{Pia) ) 



(o)' 



around the pole. Z^> and m^' are extracted with p$ = 0, 



m 



(o) 



log 



m Q + r t + Jmo + 2r t m + 1 



l + r t 



'(0) 



m 



(0) 



m 



Z<® = cosh(m( 0) ) +r f sinh(m5, 0) ). 



(38) 



(39) 

(40) 
(41) 



Imposing finite spatial momenta we determine from the speed of light and from the 
dispersion relation. The results are given by 

sinh (ml )) 



,(o) 



(0) 

rrip 



(42) 



,(0) 



cosh(m(°)) + r t sinh(m^) sinh(m^ ^ 



(o) 
nip 



m p 



(oy 



1 



m 



(0) 



(Zf-i/f)). 



(43) 
(44) 



The one-loop contributions to the quark self-energy are depicted in Fig. El whose expres- 
sion is given by 



n d 4 k 

g 2 E(p,m ) = g 2 J^—— i I J ;(k,p,mo) 



9' 



(2tt)< 

i'jo sin p B (p, m ) + vi li sinpiB^p, m ) + C(p, m ) 



(45) 



Incorporating this contribution, the inverse quark propagator up to the one-loop level is 
written as 



S 1 (p,m) = ij sin p [l - g 2 B (p,m)} + ui^2-fismpi[l - g 2 Bi(p,m)} + 



m 



+2r t sin 2 ( | j + 2r s £ sin 2 f | ) - g 2 C(p, m), 



where we redefine the quark mass as 



m = vtlq — g 2 C(p = 0, m = 0), 
C(p,m) = C(p, m) — C{p = 0, m = 0). 



(46) 



(47) 
(48) 



With this definition the inverse quark propagator satisfies the on-shell condition for the 
massless quark up to the one-loop level : S~ 1 (po = 0,Pi = 0, m = 0) = 0. For convenience 
we replace the tree level pole mass of eq. (J39~j) by 



m!^ = lo£ 



m + r t + y/m 2 + 2r t m + 1 



l + r t 



(49) 



where mi -* = at m = 0. In the following analyses we use m as if it were the "bare" quark 



mass. 



The pole mass m p is obtained from the pole of the quark propagator: 5 1 (po = im p ,pi 
0, m) = 0. We obtain 



m 



(i) 



,(0) 



Tr 



(7o + 1) 



{E(p*,m)-£(p = 0,m = 0)} 



(50) 



where p* = (po = im p °\pi = 0). It is noted that rrvp has no infrared divergence. In Fig. 0] 
we plot the dependence of mj^ for the plaquette and the Iwasaki gauge actions. We 
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observe that vanishes at m^'—0 as expected. The solid lines denote the fitting results 



of the parameterization: 



m 



(i) 



(51) 



The relative errors of this interpolation are less than 1% over the range < mjf) < 10. We 
tabulate the values of the parameters a, and 6j (i = 1, . . . , 3) in Table 

The wave function Z q is defined as the residue of the quark propagator S(p,m). The 
one-loop contribution on the lattice is given by 



= {sinh(mj, 0) ) + r t cosh(mj, 0) ) W 1} - Tr 



(7o + 1) / .x 9Siatt ( * v 



(52) 



The infrared divergence in Z)p is regularized by introducing the fictitious gluon mass A. 
We extract the divergent term in Z^ by subtracting from Is an analytically integrable 
expression is which has the same infrared behavior as 1%. As a candidate of is we take 

1 1 



Uk,p,m^) = C F Z^e(A 2 -k 2 )z lo 



with a cut-off A (< 7r). The integration is easily performed jiol] 



(53) 



-Tr (2vr) 4 ic 
v(o)_^ 

9 167T 2 



-(fc,p*,m(°); 





A 2 


-2 log 


A 2 



3 A 4 



9 A 



4 (o) 4 



2 mf 



A 2 + Amp 



(oy 



+ ^(A 2 + 4mf 2 )l-61og 



A + v // A 2 + 4mJ, 0)2 



2m P 0) 



(54) 



It is noted that the divergent term is the same as the Wilson case in Ref. jiol ]. 

The finite renormalization factor from the lattice regularization scheme to the continuum 

_ i 

Naive Dimensional Regularization (NDR) scheme (^Wt = Z q 2 ^iatt) is determined by 

Z n 



Zq = (4t t )"[i-5 2 A«K 0) : 



^g,cont 
^9,latt 



(55) 



up to the one-loop level. Here the continuum wave function renormalization factor is given 
by 



y -i 2 



(56) 



10 



with 



S C ont(P, mi 0) ; 



d u k 



1 



1 



^7o 



(57) 



where D = 4 — 2e (e > 0) in the NDR scheme. The momentum assignment is depicted in 
Fig. El(a). After some algebra we obtain 



<9£ 



cont 



(p* , m 



(0)A 

p > 



16tt 2 



-- - log 



(0)' 



2 log 



A 2 



(58) 



where 1/e = 1/e — 7 + ln(47r) and /i is the renormalization scale. In the MS scheme, the pole 
term 1/e should be eliminated. From eqs.(|52|) and ()56|) the finite renormalization constant 



AW(m™) is expressed as 



<9£ 



cont 



idp 1 



(p*,mf) 



sinh(m^) + r t cosh(m^ 



,(0) 



7 (0) 



Tr 



(70 + 1), .^Si att „ 

■«)— \p , mj 



4 



<9p 



(59) 

A=0 



Comparing eqs. (|54|) and (|58|) we find that the infrared divergence for A — > and the mass 
singularities at — ► are exactly canceled out, which assures that A^(m^ ^) is finite even 
in the massless limit. Figure shows the dependence of AW(m^) for the plaquette 
and the Iwasaki gauge actions. We parameterize as 

Eti«*K o) r 



A (1) = A (1) im 

q q \ p 



(0) 



0) + 



i + EtiM4 0) } 



(60) 



where the values of A I 1 ) (m. 



(0) 

p 



„ ..,,.„ 0) except for the DBW2 gauge action are taken from Ref. 
The fitting results are drawn in Fig. El by solid lines, whose relative errors are at most 1% 
over the range < rrvp < 10. The values of the parameters a,, b{ (i = 1,...,4) and 
Af)(mf = 0) are listed in Table HD 

The parameter v is determined by adjusting the speed of light in S~ 1 (p,m). Comparing 
the coefficients of 70 and 7, in the numerator we obtain 

sinh(m p ) [1 — g 2 B (p*, m)] 



v 



m r 



[1 - g 2 B l (p*,m)} ' 



(61) 



The one-loop contribution is given by 



'cosh(m. 



(o)> 



(o) 
rrip 



sinh(m^ ( 



(o)> 



m 



(i) 



+ uW{B i (p*,m)-B (p*,m)}, 



(62) 
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where Bi and B have no infrared divergence. The quark mass dependences of u^/u^ for 
the plaquette and the Iwasaki gauge actions are shown in Fig. |H1 As expected vanishes 
at m^ = for both cases. The solid lines depict the results of the interpolation: 



(63) 



The relative errors of this interpolation are less than a few % over the range < mi ' < 10 



The values of the parameters a, and bi (i = 1, . . . , 5) are collected in Table IITTI 

The parameter r s is determined from S~ 1 (p,m) such that the correct dispersion relation 
is reproduced: 



(64) 



This condition yields 
1 



(o) 
nip 



' (0) W, 



m) 



!/«} 



m 



(i) 

E_ r (o) 

(or * 



+ Tr 



;i + 7o) dY, 
2 dpi 



{p*,m) 



The infrared divergence of the last term can be extracted by using in eq. (|53|) : 

d 4 k dh 



n (2tt) 4 dp) 



-(k,p*, m. 



(0)N 



m P 0) 16vr 2 



log 



A 2 



A 2 



A 2 



1 A 4 



(0)' 

rrip 



2 ml 0) 



- 2- 



(o) 
m p 



-^V A2 + 4 ^ 0)2 



^(A 2 + 4mf 2 )l-21og 



ml 



A + yA 2 + Anip 



2m P 0) 



(65) 



(66) 



We find that the infrared divergence and the mass singularity in the last term are exactly 
canceled out by those in Z^/nvp. We show the m^ dependence of r^/r^ ' for the pla- 
quette and the Iwasaki gauge actions in Fig. [71 where r)P / r^ ' becomes close to zero as m^ 
vanishes. This is an expected behavior because the deviation of r s from r t stems from the 
power corrections of m^a. The quark mass dependence of r^/r^ is well described by the 
interpolation with the relative errors of less than a few % over the range < m(°> < 10, 



^ 0) l + £f=i6i{4 0) } 



(67) 



We give the values of the parameters a,i and bi (i = 1, . . . , 5) in Table KV\ 
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V. DETERMINATION OF c E AND c B UP TO THE ONE-LOOP LEVEL 



We employ the on-shell quark-quark scattering amplitude to determine eg and c#. At 
the tree level the parameters u, r s , ce and c# are adjusted to reproduce the continuum form 
of the scattering amplitude at the on-shell point removing the rriqa corrections^], 

T = -g 2 (T A ) 2 u(p')-f ll u(p)D llu (p - p)u(q'huu(q) 
-9 2 (T A ) 2 u(q')^ ll u(p)D flu (p - q')u{p')j u u(q) 

+0((p l a) 2 ,(q l a) 2 ,(p' i a) 2 ,(q' i a) 2 ), (68) 

where the momentum assignment is depicted in Fig.|H]and D^ u denotes the gluon propagator. 
At the tree level the quark-quark-gluon vertex is written as 



u(p')A$\p,p')u(p)) iatt = Zf W)i lQ u{p)) COTA + 0(( Pl a) 2 , {p\a)% (69) 
u{p')kf(p,p')u{p)\ att = Zf (u(p')*lku(p)) cont + 0((p t a) 2 , {p' t a)% (70) 



for 

A (o), f\ ■ (Po+Po\ , ■ (Po+Po\ 
Ao \P,P ) = «7o cos I — - — I + r t sin I — - — I 

(o) / _ / \ 

+ f|_ cos ^£2_P»j^^ sin(! ,,_ ! ,; ) , (71) 
A«»W) = iv m JtCOS (Et±I^ +r (». sin jEi±iiJ 

(0) / / \ 

+ — cos I — - — I a k0 sm(p - P ) 

(o) / _ I \ 

cos E MPi ~ Pi), (72) 

\ / l^k 

where the spinor on the lattice is given by 



u{p)= | | +0(( Pt a) 2 ), (73) 

with N(p) = (— i)sin(po) + mo + r t (l — cos(p )) + r s Z)i(l — C0S (Pi))- The 0(a) improvement 
condition yields 



,(») = ( 74 ) 
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.(0) 



.(0) 



JO) 
C B 



cosh(m| ) - ) ) + r t smh(m^] 



sinh(mj, ( 



(0)> 



m 



(0) 



(0)' 



r t z/ 

r (0) 



(0) 



(75) 

(76) 
(77) 



It should be noted that the values of z/(°) and r|°) are exactly the same as those determined 
from the quark propagator. 

Let us turn to the one-loop calculation. Recently the authors have shown the validity 
of the conventional perturbative method to determine the clover coefficient csw U P to the 
one-loop level in the massless case from the on-shell quark-quark scattering amplitude 3] . 
We extend this calculation to the massive case. According to Ref. jjjj], it is sufficient for us 
to improve each on-shell quark-quark-gluon vertex individually. To determine the one-loop 
coefficients c|; and c B we need six types of diagrams shown in Fig. El We first consider 
to calculate Cg . Without the space-time symmetry the general form of the off-shell vertex 
function at the one-loop level is written as 



Aj; 1} (p,g,m) 



i=a,...J 

E 

i=a,...J 



d 4 k 



where A k (p, q, m) = (p, q, m 



* (2tt) 4 " 

7**i + lk {j>F k + o F k } + {iFt + i,F k } lk 

+MF k + M F 7 k + i^F k 

+ (p k + q k ) [#* + i>H k 2 + 4H k + HH\ 

+(p k - q k ) [G\ + iG\ + 4G\ + 0G k ] + 0(a 2 ), 

g 2 A k 1) (p,q,m) + 0(g A ) and 



(78) 



^ , PoPieti 
a=0 
3 

a=0 

Polo, 

qoio- 



(79) 

(80) 
(81) 
(82) 



The coefficients F„- 



1 8), G; 



1,2,3,4) and iff (i = 1,2,3,4) are functions of 



p 2 , q 2 , p-q and m. From the charge conjugation symmetry they have to satisfy the following 
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condition: 

F k = Ft (83) 

F k = Ft (84) 

F 7 k = Ft (85) 

H k = H k , (86) 

G k = G\ = 0, (87) 

G\ = -G\. (88) 

Sandwiching (p, q, m) by the on-shell quark states u(p) and u(q), which satisfy ^u(p) = 
im p u{p) and u(q)^ = im p u(q), the matrix element is reduced to 

= u{q)lku(p) {f? + tm p (F k + F k ) - m 2 p F k ) 
+u(q)jk r Tou(p)(po ~ <?o) {F k + im p Ff} 
+ (p k + qk)u(q)u(p) {H k + im p (H k + H k ) - m 2 p H k } 

HVk - q k )u(q)u(p) {G k + im p (G k + G k ) - m 2 p G k ) + 0(a 2 ), (89) 

where we use F k = F k and F k = F%. (Note that we can replace m p with rrv^ in the 1-loop 
diagrams.) The first term in the right hand side contributes to the renormalization factor 
of the quark-quark-gluon vertex, which is equal to at the tree level. From eqs. ()87|) and 
(JHEJ we find that the last term of eq. ()89j) vanishes: this term is not allowed from the charge 
conjugation symmetry. It is also possible to numerically check G k +im p (G k +G k )—m 2 G k = 0. 

The contribution of the second term is 0(a 2 ). This can be shown as follows. For simplicity 
we first consider the case of Aqcd <C rriQ <d _1 . The difference of po and go is expressed as 



(po - q )a = [^Jml + p 2 - ^Jm 2 p + q 2 ) 



= (^)a + o(&a,£a). (90) 

Since the terms jk$o > 4 oik, 4lk]/>o and ^o7fc^ represent the violation of the Lorentz symmetry 
due to the finite rriQa corrections, their coefficients should vanish at the massless limit, 
namely F k , F k cx m p a as their leading contributions. Hence the combination of F k , F k and 
(Po ~ qo) a results in 0(a 2 ). This is retained even in the case of Aq CD -C ttiq ~ a -1 . 
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The relevant term for the determination of eg is the third one, which can be extracted 
by setting p = p* = (p = im p ,pi = 0) and q — q* = (q — im p , q^ = 0) in eq. (J78*j) : 



Hi + %m v {H k 2 + Hi) - m%H\ 



p=p ,q=q 



i T r 



() + Al A M( p *,g*,m)( 74 + i) 



dp k dq k 



-Tr 



[(- 


--) 


[dpi 


dqi J 



A^{p*,q*,m)(j 4 + l)j ilk 



i^k 



(91) 



where we have used the fact that F k , G and H k are functions of p , q and p ■ q, so that 



dFf 






dpi 


p=p*,q=q* 


dqi 


dHf 




dHf 


dpi 


p=p*,q=q* 


dqi 


dGf 




dG k 


dpi 


p=p*,q=q* 


dqi 



p=p*,q=q* 



p=p ,q=q 



p=p* ,q=q* 







(92) 
(93) 
(94) 



with j = 1, . . . , 8, I = 1, 2, 3, 4 and % = 1, 2, 3. 

We should remark that the third term in eq.(j89j) contains both the lattice artifact of 
0(p k a,q k a) and the physical contribution of 0(p k /m,q k /m). The parameter eg is deter- 
mined to eliminate the lattice artifacts of 0(p k a,q k a): 



_ r (i) 



H\ + im v {H\ + H k z ) - mlH k A ^ 
1 p\ i si p i ip=ptq=qt 

H\ + im p (H* + Hi) - mX n °" :i 



p=p ,q=q 



(95) 



where we take account of the tree level expression for the quark-quark-gluon vertex in eq. (|7U|) 
and eq.(3.51) in Ref. Q. 

We first show the calculation of eq. ()9H) in the continuum theory. The contributions of 
Figs. M (a) and (b) are expressed as 

d D k 



(2n 



icltik, P,q,m) 



with 



r(a) 



'fc,cont(^?P5 9) m ) 



1 



1 



2NJ , W + ^) +m (' 



(o) 



ilk 



1 1 



(96) 



(97) 
(98) 



x [S a p(-Pk + Qk - 2k k ) + 8 k p{~Pa + q a + K) + 5 ka {2pi3 - 2qp + M] 
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where we have replaced m p with m^p in the 1-loop diagrams. Note that Figs. El (c), (d), (e) 
and (f) do not exist in the continuum. Applying the formula of eq. (j91j) we obtain 



1 



2NJ 16vr 2 



m 



(0)' 



1 1 



16vr 2 



(0) 



log 



m 



A 2 



+ 3 



(99) 
(100) 



Here it should be remarked that we find the same results for the time component of the 
vertex function because of the space-time symmetry in the continuum theory. 

To investigate the infrared behavior of the lattice integrand I^j att (k,p, q, m) we expand 
it in terms of k. The following terms possibly yield logarithmic divergences: 



(b) 



(c) 



2N t 
Nc 
2 



2mf(c^-rf)J a (k^\X), 

-mf\cP + 2r^)J b (k,m^,X) 

+4Zf {j b (k, mf\ A) + m^J c (k, mf , A)} 

V, 
4 



(S)cPj d (k,X), 



where 



J a (k,m^,X) 
J b (k, m p,\) 
J c (k,mf\X) 



1 



1 



k 2 + X 2 (k 2 + 2im ( ^k A f 
1 iki 

(P + A 2 ) 2 p _ 2^%' 
1 fc? 



(k 2 + A 2 ) 2 (k 2 
1 



2im^kA 2 ' 



(101) 

(102) 
(103) 

(104) 
(105) 
(106) 
(107) 



(k 2 + A 2 ) 2 

with no sum for the index i. Figures El (d), (e) and (f) have no infrared divergence as long 
as m ^ 0. The coefficients of the logarithmic divergence for Jj (i = a, b, c, d) are obtained 
by performing the integration with the cutoff A: 



* d A k 

,j2^y 



6(A 2 -k 2 )J a (k,m p °\X) 



d A k 



9{K 2 -k 2 )J b (k,m^,X) 



16tt 2 
-log 
1 



,(0)' 



^log 



A 2 



A 2 



mp 

A + VA 2 + Amp 



(0) 2 



2m 



(0) 



167T 2 m f 





A 2 




A 2 


-2 log 



1 1 A 2 



(108) 
(109) 
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1 A 
+- 



4 (o) 
rrip 



?v A 2 + 4m( 0) ' + lof 



A+ VA 2 + 4m. 



2m p 0) 



d A k 



tt (2tt) 



0(A*-^)J c (A;,mf ,A) 



2m 



(0) 



6(A 2 -k 2 )J b (k,m^,\), 



d A k 



-■n (2?r) 



^(A 2 -P)J d (A;,A) 



1 




A 2 






log 




- 1 


16tt 2 


A 2 





(110) 

(111) 



From these results we find the coefficients of the infrared divergence in Figs. |H1 (a), (b) and 
(c): 



(b) 

(c) 



2N ( 

Nc 
' 2 

Nc 
' 4 



(4 0) 



(o 

s L 



trip 



-3)4 0) L 



(112) 
(113) 
(114) 



where 



16tt 2 



In 



A 2 



A 2 



(115) 



Taking the summation the total contribution is 



2N, 



,(°) 



zm_ 

(0) 



L. 



(116) 

we are left with 



If the tree level values are properly tuned as c 
-(-N C /2)(Z^ /m^)L, which is exactly the same as the infrared divergence in the con- 
tinuum theory with the correct normalization factor. 

Figure [TUl shows the quark mass dependences of Cg/ c§' for the plaquette and the Iwasaki 
gauge actions. We find relatively modest quark mass dependences for both cases. The solid 
lines denote the fitting results of the interpolation: 



c (1) 

C B 

c (0) 



c (1) 

C B 

c (0) 

C B 



+ 



,(°)- 



(117) 



where the values of the parameters Oj and h (i 



with c$ / c$ at m p °) = taken from Ref. |3j. The relative errors of this interpolation are 
less than a few % over the range < < 10. 



5) are given in Table ITTT1 together 
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We now turn to the calculation of c% . The general form of the off-shell vertex function 
for the time component at the one-loop level is written as 



A { o\p,q,m) 



12 A o l \p>(i> m ) 

i=a,...,f 

12 / TTT-U 1 ^ (k,p,q,m) 
loF? + 70^2 + ^70^3° + MF, 

+( Po + go) [^ 1 + ^ 2 ° + ^3° + 
+(Po-<?o) \G*i+i>Gl + iGl- 



0(a 2 ), 



(118) 



where we define A (p, q, m) = Ao°^(p, q, m) + g 2 AQ L) (p, q, m) + 0(g 4 ). The coefficients Ff, G° 
and H® (i = 1, 2, 3, 4) are functions of p 2 , q 2 , p ■ q and m. As in the case of the charge 
conjugation symmetry provides the coefficients with the following constraint: 



Fl 



F° 

#3> 



— ^4 



G° 2 



(119) 
(120) 
(121) 
(122) 



Sandwiching Aq 1 (p, q,m) by the on-shell quark states as before, the matrix element is re- 
duced to 



u{q)A { o\p,q,m)u(p) 
u{q) lQ u(p) {F? + im®(F° + F 



m 



(0) 



+ (Po + q*)u{q)u{p) {fl? + imf\H° + H° 3 ) - m^H*} 

+ (Po ~ qo)u(q)u(p) {G? + imf (G§ + G§) - mf 2 G°} + 0(a 2 ). (123) 

Here we replace m p with The renormalization factor is determined from the combi- 

nation of F® + imf\F% + F 3 °) — m^ 2 F°, which should be the same as F^ + im^(F£ + 
in eq. (j89jl . The last term of ecj. (J123J) vanishes from ecjs. (J121J) and (jl22jl as a 
consequence of the charge conjugation symmetry. We can also check G® + irrv^'(G% + Gf) — 



rn 



(fi) 2 QQ = o numerically 



The coefficient eg is determined to remove the 0(a) contribution of the second term in 
the right hand side, where the physical contribution of O(p /m, qo/m) is also included. The 
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second term is extracted by setting p = p* and q = q* in eq. 



as 



p=p*,q=q* 



\{-- 


h— ) 


_{dpi 


dqi J 



• (0) 

imp 



1 iTr 



A ( a\p*,q*,m) 



-Tr 





--) 


.1% 


dqi J 



Ao 1) (p*,g*,m)7 i 7 4 



-zmi 0) -Tr 
p 4 



d 2 



Ao 1) (P*,?*,^)(74 + l)7i7i 



(124) 



dpidqj 

where we again have used the fact that F°, G° and H° are functions of p 2 , q 2 and p ■ q. We 
determine the parameter ce to remove the 0(pka,qka) contributions: 



^ - ,y(l) 



p=p* ,q=q* 
2 



-Zf H[ + imp) (H° 2 + H° 3 ) - mf X° 



cont 

p=p*,q=q* 



(125) 



where we take account of the tree level expression for the quark-quark-gluon vertex given in 
eq.fU) and eq.(3.50) in Ref. Q. 

The infrared behavior of the integrand Iq (k,p,q,m) is examined by expanding it in 
terms of k. The logarithmic divergences are attributed to 



2NJ {zy(o) 2 + r (°) s inh(m^ 0) )} 2 
(°) 2 {4i/°) sinh(mf) - 3^°) Sinh( ^ 0)) - m<V > a }Zj 0) 



m 



+2z/°){m( ( V ) - sinh(mf)}{z/(°) 2 + 2rf sinh(m( 0) )}ZW 



+{- 
+{- 



sinh(m 



(0)> 



(0) 

trip 



+ m (o) i ,(o) 2 }co S h(m(°))Z( ) : 



sinh(m^°)) 2 



m 



(0) 



m£V°) }r t sinh(mf ^ 



(6) 



+2mf^cf sinh(mf )Zf 2 ] , 
jVA J 6 (fc,m(°),A) 



(126) 



2 / ^(o) 2 _|_ r (°) sinh(mp°)) 



[-24 0) sinh(m(°)){m(°) r (o) + z ,(o)} 



. C (0) Z ( ) ^(0)^0) +sinh(m (o )) | 
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(0) 



nip J 



-V 0) {sinh(m( 0) ) - z/ (0) m( 0) } {sinh(m( 0) ) + r t cosh(m p 0) )} 



+ - 



iV, 



J c (fc,m(°),A) 



2 7 ^(o) 2 _|_ r (°) s inh(m 



,(o)> 



f 4 0) sinh(m(°))Z( ) 



-2Z/(°)Z(°) {sinh(mf ) + 3mJV )} 



^) (-3)4°>J d (M)- 



(127) 
(128) 



We find no infrared divergence for Figs. El (d), (e), (f) as long as ^ 0. The momentum 
integration with the cutoff A yields the following logarithmic divergences: 



(b) 



(c) 



2N c J 2mf {^°) 2 + rf sinh(mj 0) )} 2 



x 



z/(°) 2 {4z/W sinh(m(°))_3z/°) 



sinh(m| ) - ) ) z 



(o) 



+2z/°W ( V 0) - sinh(m(°))}{z/(°) 2 + 2rf sinhrf)}^ ) 



. g^^j + ^o)^) 2 } cosh(m (°))^ 



(0) 



+{ _sinh(mp_ _ m (o )z/ (o) 2 }rtSinh(m (o ))z (o)^ 
m p 

+2m(°) z/ (0)4 ) sinh ( m (0) )Z (0) ; 



(129) 



-1) 



(0)v 



2 7 2mp°- > z/(°) 2 + ri ' 1 sinh(mj 
|4 0) Z(°){ i ,(0) m (0) +sinh(m j)) ) j 

+„(o) f_ 3 sinh ^i 0) ) + ^ojj ^(0) + r (o) m (o) j 



-24 0) sinh(mf ) {mfrf W 0) } 



(°) 

--z/W {sinh(mf) 



Am P 0) z/(°) 2 + rrsinh(mH 
-Iv^Zf {sinh(m(°)) + 3mJV°>}' 

-f) (- 3 )4 0) ^ 



{sinh(m( 0) ) + r t cosh(mj, 0) )} 



-8mf 4 0) sinh(mf )Zf 



(130) 
(131) 



Once we demand the tree level conditions, 



r (o) 



(132) 
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+ rf sinh(mf) = Zf\ (133) 



7/ v 



, rn sinh(m^) 
u (0) = Kp_f_ } (134) 

rrip 



the above expressions are reduced to be 



(a) : (135) 

(b) : (-^) * L {-2mp)c^ + 3c£><> - 2z/°>} (136) 

(c) : (-^ (-3)4 0) L. (137) 



4 

Finally, with the aid of another tree level condition 



'b = To) > (138) 



the total contribution is found to be 

which is the same as that for the space component in eq. (|llfi|) . 

We again stress that the infrared divergences originating from Figs. El (a), (b), (c) contain 
both the lattice artifacts and the physical contributions. The former exactly cancels out if 
and only if the four parameters r^°\ c$ and Cg are properly tuned as denoted in 
eqs.(J22J), (|75|) . (|75|) and (|77|) . This is another evidence that the tree level improvement is 
correctly implemented in Ref. 

In Fig. [TT]we plot Cg / ' c$ as a function of for the plaquette and the Iwasaki gauge 
actions. The fitting results of the interpolation 



r (o) ~ Jo) 



+ 7 %fd^ (140) 



4 0) =o 



1 + ELiWT 



are also shown by the solid lines. We find relatively modest quark mass dependences similar 
to the cb case. The relative errors of this interpolation are less than a few % over the range 
< mf < 10. Table EH summarizes the values of the parameters aj and bi (i = 1, . . . , 5) 
and c$ / Cjp at = taken from Ref. jjj. 
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VI. MEAN FIELD IMPROVEMENT 



In this section we rearrange the 1-loop results in the previous sections, using the mean- 
field improvement. We first replace the link variable U nit by u(U nii /u) = uU n ^, where u is 
the average of the link variable u = (U n ^) in some gauge fixing, or u = P 1 / 4 with P is the 
average of the plaquette. In this paper, we adopt the latter definition: 

u=l-g 2 ^-T MF . (141) 



A detailed description on the derivation of Tmf is given in Sec. Ill of Ref . pjj • 
This replacement leads to the following dispersion relation. 

usmh(rh^) = m + r t (l — wcosh(m^)) + 3r s (m p °')(l — u) 

= m + r t u{l - cosh(mj, 0) )) + (1 - w)3(r s (mj, 0) ) - 1) (142) 

where m = mo + r t (l — u) + 3(1 — u) and fh p ^ is the tree level pole mass in the mean-field 
improvement. Using the relation 

sinh(m^°- ) ) + r t cosh(m p ^) = m + r t , (143) 

we have 

(1 + u— l)(sinh(m^°- ) ) + r t cosh(m£°))) = sinh(m^) + r t cosh(m^°- ) ) 

+ r t (u - 1) + (1 - u)3(r a (m< 0) ) - 1). (144) 

This leads to the relation: 

(o) ~ (o) i / ^ sinh(m(°)) + r 4 (cosh(m(°)) - 1) + 3(r s (mf ) - 1) 

cosh [nip ) + r t smn(mp J 
= m p 0) + g 2 Am p , (145) 

where 

C F sinh(m(°))+r t (cosh(m(°))-l) + 3(r s (m( ))-l) 

Am P = 2~ ■ w~(oh • ( 146 ) 

z cosh(mp ) + r t smh(m,p ) 

Using this, the pole mass is rewritten as 

m p = mf + g 2 m p V = fhf> + g 2 m p ^ (147) 
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where mj, 1 ' = mj, 1 ' + Am p . 

With the use of fhf£> we apply the mean field improvement to Z q , u, r s , ce and cb- 



^,latt - z 9l iattl m p J« 1 + # — (o) r- fi- — Jmf + — (0) (0) Am, , (148 J 



g,latt g,latt 

1/ = ^\mf) + jV^mf) + g 2 ^Am p , (149) 

amp 

r s = rf(mf ) + ^ 2 r«(mf) + &Am p , (150) 

amp 

1 / ^ q ^2 Q.(0) \ 

(0) 1 / -, . 2 C E 2 6 n rp 9 OC E A I /-.r^ 

Ci7 = c * V + ^ 1°) ~ 9 2 Cf mf + ^y^lsy Amp ' ( 151) 

c B = c B >-^l + g -^-g -C F T MF + -^^Am p J. (152) 

One then finally replaces u = P 1 / 4 with the one measured by Monte Carlo simulation. 

The mean-field improved MS coupling g^(fi) at the scale \i is obtained from the lattice 
bare coupling g^ with the use of the following relation: 

IP 22 / 4 \ 

= ig + *' + Cp + 16^ l0g(/ia) + ^ - 48^ l0g(/ia) J • (153) 

For the improved gauge action one may use an alternative formula [l^ 
1 c P + 8c x Pl + 16c 2 P2 + 8c 3 P3 
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+d g + (co ■ c p + 8ci • c R i + 16c 2 • c R2 + 8c 3 • c m ) + ——: log(yua) 

l07T z 

+]V /( d /-^ lo sW)' (154) 

where 



p = 


\™ pl = 1 - 


-c p gl + 0(g^), 


(155) 


Rl = 


-TrU rtg = 1 


-cm g% + 0(go), 


(156) 


R2 = 


-TrU chr = 1 


- c R2 g 1 + O(g Q i ), 


(157) 


R3 = 


-TrU plg = 1 


- c R3 gl + 0(go), 


(158) 



and the measured values are employed for P, Rl, R2 and R3. T 
crz for various gauge actions are listed in Table XVI of Ref. jll 
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re values of c„, cri, cr 2 and 



VII. CONCLUSION 



In this paper we determine the 0(a) improvement coefficients, is, r s , cb and eg in the 
relativistic heavy quark action up to the one-loop order for the various improved gauge 
actions. As byproducts we also calculate the quark wave function Z q and the pole mass 
m p . While v, r s , Z q and m p are determined from the quark propagator, we use the on-shell 
quark-quark scattering amplitude for eg and eg. The mqa dependences are examined by 
making the perturbative calculations done in the rriQa dependent way: As for the results of 
i/CO/iA ), rf'/rf), Z^/Z^ and rrvp we observe the strong mqa dependence for rriqa-^l, 
while the dependence becomes much milder beyond mqa ~ 1. On the other hand, eg / Cg 
and c^/c^ show relatively mild mqa dependences for < mqa < 10. Employing the 
conventional perturbative method with the fictitious gluon mass to regularize the infrared 
divergence we show that the parameters is, r s , cb and ce in the action are determined free 
from the infrared divergences. This is achieved if and only if the tree level values for is, r s , 
cb and ce are properly adjusted as presented in Ref. For later convenience we give a 
detailed description about how to apply the mean field improvement to our results. We are 
now trying a numerical test of this formulation with the mean field improved parameters 
employing the heavy-heavy and heavy-light meson systems. 
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TABLE I: Values of parameters a% and bi (i = 1, 2, 3) in the interpolation of m p with ea.(|51|) for 
the various gauge actions. 



gauge action a\ a<i a% b\ b<i b^ 



plaquette 


0.44498 1.1694 0.20262 5.1026 3.4713 1.3421 


Iwasaki 


0.33282 1.2533 0.24384 7.7008 7.4235 2.4353 


Symanzik 


0.39766 1.2161 0.22280 5.9914 4.6478 1.6915 


Iwasaki' 


0.33819 1.2549 0.23934 7.5388 7.1244 2.3482 


Wilson 


0.32642 1.2626 0.23357 7.9398 7.8574 2.5026 


DBW2 


0.25898 1.1989 0.24076 10.889 14.868 4.1077 



27 



to 

GO 



TABLE II: Values of parameters o« and bi (i = 1, ... ,4) in the interpolation of A^ with eq. (|60j) for the various gauge actions. The values 



of A { q ] at m { p ] 



are taken from Ref. 



except for the DBW2 gauge action. 



gauge action A^(m?=0) a, 



a-2 



as 



a i 



n 



plaquette -0.07773 0.18777 3.1560 -0.15124 0.090311 15.929 11.134 0.83734 0.44445 

Iwasaki -0.01478 0.088032 4.5405 -0.59509 0.12110 48.391 27.665 -1.7715 0.83723 

Symanzik -0.05044 0.13224 4.0242 -0.51558 0.20496 25.697 17.632 -0.67533 1.3005 

Iwasaki' -0.01739 0.093314 2.8033 -0.37819 0.078356 28.654 16.445 -1.1042 0.53359 

Wilson -0.01068 0.082214 3.4543 -0.40922 0.085302 39.349 21.267 -1.0475 0.56718 

DBW2 +0.02029 0.036401 0.33748 0.19697 0.027077 9.0916 6.0990 3.1026 0.11858 



TABLE III: Values of parameters a« and bi (i = 1, . . . , 5) in the interpolation of with ea.(|63|) for the various gauge actions. 



gauge action 


a\ a 2 




04 CI5 


h 6 2 h 6 4 


65 


plaquette 


0.0013787 0.23077 


2.3586 


0.41097 0.21571 


21.776 23.401 16.822 1.0718 


2.0084 


Iwasaki 


0.010681 0.10966 


1.4455 


0.29523 0.14143 


17.745 29.424 21.954 1.3323 


2.7265 


Symanzik 


0.0063804 0.18587 


1.9939 


0.37153 0.18560 


20.428 26.116 18.407 1.2431 


2.2404 


Iwasaki' 


0.010248 0.11663 


1.4875 


0.30469 0.14531 


17.882 28.950 21.577 1.3199 


2.6761 


Wilson 


0.010284 0.10466 


1.3679 


0.29382 0.13690 


16.951 28.693 22.473 1.2577 


2.7941 


DBW2 


0.011730 0.021615 0.77692 0.18009 0.085451 


13.817 30.339 31.211 0.83789 3.8899 



TABLE IV: Values of parameters and bi (i = 1, . . . , 5) in the interpolation of 7s with eq,(|67[) for the various gauge actions. 



gauge action 


ai 


«2 a 3 


a i 


o 5 bi b 2 b 3 bi 


65 


plaquette 


0.048659 


0.23362 1.6049 


0.15709 


0.068748 6.9070 18.195 11.853 0.74864 


0.92944 


Iwasaki 


0.0015650 


-0.14504 1.3871 


0.012363 


0.044884 12.493 27.761 31.031 -2.4872 


2.0163 


Symanzik 


0.024680 


0.14420 0.92264 0.11171 


0.044813 6.8434 12.778 12.079 0.010953 0.97640 


Iwasaki' 


0.012139 


-0.18772 2.3352 


0.037267 


0.080842 20.355 42.341 48.503 -3.5202 


3.1917 


Wilson 


0.015805 


-0.14494 1.9526 


-0.034214 0.066725 17.479 36.192 41.823 -3.7429 


2.5913 


DBW2 


-0.038636 


-0.21232 0.44240 


-0.11604 


0.015381 5.3026 20.770 22.879 -9.2470 


1.6088 



TABLE V: Values of parameters a« and bi (i = 1, . . . , 5) in the interpolation of c B with eq, (|117j) for the various gauge actions. The values 
of Cg^/c^ at = are taken from Ref. [J. 



gauge action 




= 0) oi 


a-2 a 3 


a. i 


a 5 


h b 2 


6 3 b 4 b 5 


plaquette 


0.26858825 


0.11031 


-0.45775 1.4585 


3.5629 


0.19832 


-1.2454 6.9947 


72.053 7.4057 6.4400 


Iwasaki 


0.11300591 


0.031328 


-0.25330 -0.97944 4.9518 


-0.70873 


11.329 -16.089 275.53 27.127 13.769 


Symanzik 


0.19624449 


0.070876 


-0.35318 0.75783 


2.6525 


-0.034442 


-2.6878 6.1987 


70.755 8.3789 5.2277 


Iwasaki' 


0.12036501 


0.025289 


-0.18878 0.14872 


1.3433 


-0.17271 


-0.73798 9.9278 


71.365 8.0592 3.6373 


Wilson 


0.10983411 


0.021005 


-0.19478 0.16783 


0.93978 


-0.15323 


-0.74965 11.415 


64.418 7.0567 2.8708 


DBW2 


0.04243181 


8.7741 


-14.110 -19.180 


1.0685 


-5.2296 


-516.57 1177.1 


76.371 472.69 37.314 



TABLE VI: Values of parameters a, and bi (i = 1, . . . , 5) in the interpolation of c E with eq. (|14fl|) for the various gauge actions. The values 
of CgVc^ at = are taken from Ref. [J. 



gauge action 




= 0) a x 


«2 a 3 


a. i 


0-5 


bi b 2 


b 3 6 4 


h 


plaquette 


0.26858825 


0.061511 


-0.52275 1.5301 


0.087055 


-0.009142 


-8.4208 27.507 


-2.5724 3.1868 


-0.15206 


Iwasaki 


0.11300591 


0.014754 


-0.16103 0.47732 


-0.093278 


-0.0044554 


-5.8133 21.985 


-1.4805 1.7759 


-0.025592 


Symanzik 


0.19624449 


0.046299 


-0.38345 1.0557 


-0.013178 


-0.0040176 


-7.9568 23.808 


-1.5669 2.4568 


-0.058986 


Iwasaki' 


0.12036501 


-1.0186 


10.033 -0.22227 


-0.57195 


0.048774 


533.27 -178.18 147.82 -22.688 0.92777 


Wilson 


0.10983411 


21.659 


-158.46 271.84 


-50.696 


1.1180 


-3836.5 8792.4 


-268.31 644.04 


-28.188 


DBW2 


0.04243181 


-0.084945 


-0.49122 -0.61820 0.13749 


-0.069162 


-1.2883 45.147 


-10.501 2.9737 


0.14473 




FIG. 1: Momentum assignment for the quark-gluon vertices. 



33 



0.3 
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FIG. 5: Aq as a function of m p . Solid lines denote the interpolation of A q with eq. (|6fl|) 
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FIG. 6: as a function of m p . Solid lines denote the interpolation of /v^ with ea. (|63j) 
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Q q' Q 

(a) (b) 

FIG. 8: Tree level diagrams for the quark-quark scattering. 




(e) (f) 
FIG. 9: Quark-gluon vertex at the one-loop level. 
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FIG. 10: eg 3 / as a function of Trip' . Solid lines denote the interpolation of c y B ' / ' c^' with 
eq.(Hnj) 
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FIG. 11: 4 1} /4 0) 

as a function of m p . Solid lines denote the interpolation of c E / ' c E with 
ea. lfTlOl 
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